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Abstract 


We compute the decay rates and H-e nergy distributions of B mesons into 
the final state H + X, where H can be any one of the 5-wave or P-wave 
charmonia, at next-to-leading order in the strong coupling. We find that 
a significant fraction of the observed J/ij), f>' and Xc must be produced 
through cc pairs in a colour octet state and should therefore be accompa¬ 
nied by more than one light hadron. At the same time we obtain stringent 
constraints on some of the long-distance parameters for colour octet pro¬ 
duction. 
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I. INTRODUCTION 


Exclusive B decays provide us with important information on the structure of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix. However, the theoretical calculation of 
absolute branching fractions is complicated by the fact that a rather detailed knowledge 
of strong interaction effects is required. The theoretical situation with regard to strong 
interaction dynamics improves as one considers more inclusive final states. At leading 
order in A QCD/^b, where is the b quark mass and A q C d the strong interaction scale, 
the totally inclusive B decay rate can be computed completely in perturbation theory. 
However, it is not necessary that the process be totally inclusive. A semi-inclusive decay 
B —> H + A" can also be treated perturbatively in part, provided the formation of the 
hadron H proceeds through a short-distance process. This is the case, if H is a char- 
monium state, because the production of a charm-quark cc pair requires energies much 
larger than A. The bound state dynamics of H then factorizes and can be parametrized. 
This statement is summarized by the factorization formula JE] 

F(B -► H + X) = Y, C(b -> cc[n} + x) ( 0 H [n ]), (1.1) 

n 

which is valid up to power corrections of order A QCD/m^c- (To this accuracy it is justified 
to treat the B meson as a free b quark.) The parameters (0 H [n\), defined in [|1|, are 
sensitive to the charmonium bound state scales of order m c v and m c v 2 , where v is the 
typical charm quark velocity in the charmonium bound state. With v 2 ~ 0.25 for J/-0 
we consider these scales to be too small to be treated perturbatively. On the other hand 
the coefficient functions C(b —> cc[n ] +z) describe the production of a cc configuration n 
at short distances and can be expanded in the strong coupling a s (/x) at a scale /i of order 
2 m c . We have expressed the decay rate in terms of several non-perturbative parameters 
(0 H [n\). The predictive power lies in the fact that these parameters are independent 
of the particular charmonium production process and hence are constrained by other 
charmonium production processes. 

Because charmonia pass as non-relativistic systems, Eq. (|1.1| ) involves an expansion 
in v 2 and the cc[n] configurations that appear in lower orders of this expansion can be 
usefully classified by 2S+1 Lj"\ where S, L and J refer to spin, orbital angular momentum 
and total angular momentum, respectively. In addition C = 1, 8 refers to a colour singlet 
or a colour octet configuration. In the present work we calculate the short-distance 
coefficients for 


n e 


/3q(1>8) 1q( 1.8) 3p(l) 

\ °1 5 °0 ) - r 0 ,1,2) 


3p(8) lp(l>8)| 


( 1 - 2 ) 


at next-to-leading order (NLO) in a s . As we discuss later, we believe that these terms in 
the velocity expansion are sufficient to reliably predict (to about 25%, barring radiative 
corrections in a s ) the decay rates into J/^, ry c , Xo,i ,2 an d the elusive 1 P\ state h c 
as a function of the long-distance parameters (0 H \n\). We fold that, given the present 
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uncertainties in the long-distance parameters, the experimentally observed branching 
fraction [Q for J/-0 and ip' can easily be accounted for at NLO. However, we find it 
difficult to account for the observed Xci and Xci branching fractions simultaneously, 
because the a s expansion of the colour singlet contribution in Xci production turns 
out to be untrustworthy at NLO. The NLO corrections typically enhance the decay 
rate by about (20-50)% in the colour octet channels and lead to bounds on the long¬ 
distance parameters, which should be useful for the phenomenology of other charmonium 
production processes. We also compute weights of the charmonium energy distribution, 
which yield additional information. The shape of the energy distribution itself, however, 
is difficult to predict, because it is distorted by the motion of the b quark in the B meson 
and the energy taken away in the hadronization of a cc state n. This distortion averages 
out in weighted sums, as long as the weights are sufficiently smooth. 

We then compare the inclusive calculation for J/ip+ X with the sum of the measured 
decay rates for J/ip + K and J/ip + K*. The comparison suggests a significant fraction 
of multi-body decays, consistent with the energy spectrum observed by CLEO [||. A 
substantial contribution from multi-body decays is also reassuring from the point of 
view of validity of the theoretical calculation. Factorization implies that a cc state n 
hadronizes into a J/ip plus light hadrons independent of the remaining decay process up 
to corrections of order Aqcd/^&,c- If n refers to a colour octet state, the conversion into 
charmonium requires the emission of at least one gluon. Although colour reconnections 
with the spectator quark in the B meson must eventually occur, we expect a charmonium 
produced through a colour octet cc state to be accompanied by more than one light 
hadron more often than for a colour singlet cc state. Since we find that a large fraction 
of the total decay rate is from colour octet intermediate states, we also expect a large 
fraction of multi-body final states. This evidence also suggests to us that the energy 
released in the B meson decay into charmonium is already large enough for an inclusive 
treatment to be applicable. 

Inclusive production of S wave charmonia has been considered in Refs. in the 
colour singlet model and at leading order (LO). In addition, the colour singlet production 
of the P -wave state Xci was computed in Ref. ||. (At LO the states Xco and Xci are 
not produced.) The colour singlet model is contained in Eq. ( |1.1| ) as the term where the 
quantum numbers of n match those of the charmonium state. For P -wave charmonia 
the colour singlet model does not coincide with the non-relativistic limit v —> 0 and is 
generally inconsistent. The authors of Ref. Q noted that the contribution from cc[' 3 5 < | 8 ' 1 ] 
is leading order in v for XcJ and that cc^Sq ] is leading order for h c . They computed 
the relevant short-distance coefficients to LO in a s . In the case of J/ip, the short- 
distance coefficients of cc[n] states with n = ' s S[ S \ ' i Pj H> are strongly enhanced as a 
consequence of the particular structure of the weak effective Lagrangian that mediates 
b quark decay. These production channels have to be taken into account although the 
corresponding long-distance matrix elements are suppressed by a factor of u 4 . The 
relevant coefficient functions were computed in Ref. [7J. again at LO in a s . Ref. || adds 
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a study of J/ip polarization effects. The only NLO calculation of charmoninm production 
in B decay is due to Bergstrom and Ernstrom |J , who computed the contribution of the 
colour singlet 3 Si intermediate cc state to J/ij) production. We repeated their calculation 
and comment on it later on. 

The paper is organized as follows: In Section [TT] we introduce notation and discuss 
the structure of important contributions to a given charmonium state. Section H pro¬ 
vides some details on the calculation related to the handling of ultraviolet and infrared 
divergences at intermediate stages. Section [TV| contains our main results. We present 
expressions for the decay rates in numerical form and a comparison with existing exper¬ 
imental data. Analytic results for the decay rates and energy distributions are collected 
in two appendices for reference. Section |V] contains our conclusions. 


II. PRELIMINARIES 


The terms of interest in the A B = 1 effective weak Hamiltonian 

G f 


H eff ~ ~j= ' 

V * q=s,d 


v; b v cq 


icfoMCM/*) + c [ „MO,w] - v^jzcMOM 


contain the ‘current-current’ operators 

Oi = [87„(1 - 7 s)c] [iVU - 75 ) 9 ] 

C 8 = [cT A 7^(1 - 75 )c] [bT A Y(1 - 75 ) 9 ] 


( 2 . 1 ) 


( 2 . 2 ) 

(2.3) 


and the QCD penguin operators (9 3 _ 6 . (See the review Ref. |]T0| for their precise defini¬ 
tion.) For the decays B —> charmonium + X it is convenient to choose a Fierz version 
of the current-current operators such that the cc pair at the weak decay vertex is either 
in a colour singlet or a colour octet state. The coefficient functions are related to the 
usual C± by 


q 1] (/0 = 2C + (/0 -c_(/0, 

(2.4) 

C [ 8 ] (/i) = C + ( y u) + 0 _(/i). 

(2.5) 


The NLO Wilson coefficients C±{/i ) have been computed in Refs. |[TT1 . |T2f| . With the 


conventions of Ref. 12 


C± (At) 


a s (M w ) 

(h ) 
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with 
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(2.7) 
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7± Pi r± 


Wl 
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and the one-loop and two-loop anomalous dimensions 

7± 0) — =t 2 (3 =F 1), 

'/!' = ±2 (-21 ±i n/ - 2/3l 0 «±) . 
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(2.8) 

(2.9) 

( 2 . 10 ) 


The quantity k± is scheme-dependent and depends in particular on the treatment of 75 . 
In the ‘naive dimensional regularization’ (NDR) scheme, k± = 0; in the’t Hooft-Veltman 
(HV) scheme, k± = =|=4. I 11 the HV scheme the current-current operators, implied by the 
convention used in Refs. ||T 0 |,|T 2 ||, are not minimally subtracted. If one computes the low 


energy matrix elements of the weak Hamiltonian in the modified minimal subtraction 
(MS) scheme, as we will do below, one has to apply an additional finite renormalization. 
This amounts to multiplying the coefficients C±(ji) by a factor of 1 — Aa s (fj) / (3 tt) , or, 
equivalently, to an additional contribution to k± in the HV scheme. No additional 
renormalization is required in the NDR scheme. At NLO the strong coupling is given by 




47T 


Po l 11 ^ 2 /-^- qcd ) 


Pi hi[ln(/i 2 /Ag C . D )] 
00 /^-qcd) 


( 2 . 11 ) 


with 

/So = 11 - | n f , Pi = 102 - y n f . (2.12) 

The NLO QCD corrections involve the one-loop virtual gluon correction to b —> 
cc[n] + q and the real gluon correction b —> cc[n} ± q + g, where the cc pair is projected 
on one of the states in (|1.2| ). The corresponding diagrams are shown in Figs. [L] and |2] 
respectively. The decay rate into a quarkonium can be written as the sum of partial 
decay rates through one of the intermediate cc states n. At next-to-leading order the 
partial decay rates take the form 


i>] = r 0 

aJL 


c i,8]f[ n \(v) (l + <5p[n]) 

(C^gMpq) ± 2C [1] C [8] g 2 [n\(r l ) + Cf 8] g 3 [n\(r))) ] ( 0 H [n ]), 


where 


G 2 F \V bc \ 2 ml 
2167r(2m c ) ’ 


(2.13) 


(2.14) 
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FIG. 1. One-loop virtual corrections to b —> ccq. Wave function renormalizations are not 
shown. 



FIG. 2. Real gluon corrections to b —> ccq. 


and r) = 4The operators 0 H [n) are defined as in Ref. [TJ]. The LO term is 
multiplied by C'^ if n is a colour singlet state and by C)|] if n is a colour octet state. We 
also used the fact that |V^ s | 2 + V c d | 2 ~ 1 to high accuracy. The functions f[n] and gi[n} 
will be given later. The LO contribution is multiplied by a correction term 8p[n] due to 


the penguin operators in (|2.1| ). Likewise, we write the quarkonium energy distribution 


as 


dT[ 


n\ 


dx 

+ 


= T n 


c [i,s\f[n\(v) (! + 5p[n]) 8(l + r/- x) 


(C'p] gi[n](r},x) + 2C {1] C [8] g 2 [n](ri,x) + Cf 8] g 3 [7i](r], x)) ] {0 H [n]) (2.15) 


where x = 2 P-pb/mf. Note that to leading order in A QCD/ m b we do not distinguish the b 
quark mass from the B meson mass. To the order in the velocity expansion considered in 
this paper, we can also identify the momentum of the quarkonium with the momentum 
P of the cc pair. (The kinematic effect of distinguishing the two is discussed in Ref. ||13|| .) 
Hence x can also be identified with 2 Eh/Mb, where Eh is the quarkonium energy in the 
B meson rest frame and Mp the B meson mass. 
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We now discuss which intermediate cc states should be taken into account for the 
production of a given quarkonium H. 

J/ip, xp'\ At leading order in the velocity expansion the spin-triplet S'-wave charmo- 
nium states are produced directly from a cc pair with the same quantum numbers, i.e. 
n = 3 S'j 1 ^. At order v 4 relative to this colour singlet contribution, a ip can materialize 
through the colour octet cc states n = 3 5®, 'S'®, 3 P®, where the subscript ‘J’ implies 
a sum over J = 0,1,2. The suppression factor v 4 follows from the counting rules for 
the multipole transitions for soft gluons that convert the state n into the ip meson [lj. 
The leading order colour singlet contribution is proportional to Cm, while the colour 
octet terms are proportional to C 2 ^. Because the weak effective Hamiltonian favours the 
production of colour octet cc pairs by a large factor 

qyCfi] « 15, (2.16) 

the colour octet contributions must be included, since their suppression by v 4 ~ 1/15 (for 
J/ip) can easily be compensated. (The numbers serve only as order of magnitude esti¬ 
mates of the relative importance of the colour singlet and the colour octet contributions.) 
According to the velocity counting rules, there is a correction of order v 2 to the colour 
singlet contribution related to the derivative operator ( 3 /S' ( ^ 1 ^) as defined in Ref. [1J. 
Because it is multiplied by the small coefficient Cm, and because we will find that the 
colour singlet contribution is indeed a small contribution to the total production cross 
section, we do not consider this additional correction in what follows. Similar derivative 
operators contribute to the colour octet channels. We do not take them into account, be¬ 
cause we do not take into account other corrections of order v 6 with the large coefficient 
Crg,. Hence, even after including the NLO correction in a s , there remains an uncertainty 
of order v 2 ~ 25% in the theoretical prediction, assuming that the long-distance matrix 
elements were accurately known. 

rj c . The same discussion applies to the spin-singlet state. The colour singlet contri¬ 
bution involves n = (S'®. At relative order v 4 , rj c can be produced through the colour 
octet cc states n = -(S'®, 3 S[ 8 \ 'P®. 

Xcj : At leading order in the velocity expansion, both n = 3 P® and n = 3 S'® 
contribute to the production of a the spin-triplet P- wave state 0 . Because the partial 
production rate through the '(S'® state is already multiplied by the large coefficient CL , 
it is not necessary to go to higher orders in the velocity expansion. Note that, because of 
the V — A structure of the weak vertex, a cc pair cannot be produced in a 3 P 0 ,2 angular 
momentum state at LO in a s . 

h c \ The same discussion as for the \cj states applies to the spin-singlet P-wave 
state. In this case we take into account n = W® and n = (S® at NLO in a s . Owing to 
the V — A structure of the weak vertex, a cc pair cannot be produced in a 4 P\ angular 
momentum state at LO in a s . 
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III. OUTLINE OF THE CALCULATION 


The Feynman diagrams shown in Figs. [1] and [2] are projected onto a colour and angular 
momentum state as specified in ( |1.2| ). The virtual corrections contain ultraviolet (UV) 
divergences, which can be absorbed into a renormalization of the operators (9 18 in the 
weak effective Hamiltonian (|2.1| ). The virtual corrections contain infrared (IR) diver¬ 
gences, which cancel against IR divergences in the real corrections. In addition, the real 
corrections contain IR divergences due to the emission of soft gluons from the c or c 
lines, which do not cancel with IR divergences in the virtual correction, if the cc pair is 
projected on a P -wave state. These IR divergences can be factorized and absorbed into 
a renormalization of the non-perturbative matrix elements ( 0 H [n ]). In the following we 
provide some details on the UV and IR regularization, which are specific of the present 
calculation. More details on the strategy of a next-to-leading order calculation can be 
found in Ref. [Ill], which deals with quarkonium decay and total quarkonium production 
cross sections in fixed-target collisions. 


A. UV regularization and the treatment of 75 


The UV divergences are regulated dimensionally and the IR divergences are regulated 
with a gluon mass. The UV divergences in the diagrams of Fig. [I] cancel against the 
UV divergences in diagrams (not shown in the figure) with the insertion of the 1-loop 
counterterm for We combine the diagram with its counterterm diagram before 

projecting on a particular cc state n, and before taking the 2 -particle phase space integral. 
This has the advantages that it avoids extending the projection to d dimensions and that 
the phase space integral can also be done in four dimensions. 

The finite part of the virtual gluon correction depends on the prescription for handling 
75 in d dimensions. This has to be chosen consistently with the one used to define the 
operators in Ref. [|T3- The prescription consists of a definition of 75 and its anti¬ 
commutation property, together with a choice of ‘evanescent operators’. The evanescent 
operators are implicitly defined by specifying the order e (where d = 4 — 2e) terms of the 
following products of Dirac matrices: 


7 p 7 aT Al ® 7 p 7 q= (16 + AX R e) Tp ® + E x (3.1) 

Tp7p7a ® 7c*7p r p = (4 + 4 Y R e) Tp ® + E y (3.2) 

Tp ® 7p7«Tp7 Q 7 P = (4 + 4 Z R e) + E z . (3.3) 

(Here we defined Tp = 7 p(l — 75 )-) The renormalization conventions of Refs. |p~2| , |T5f 

correspond to: 

NDR scheme : Andr = — 1 Vndr = ^ndr = —2, (3-4) 

HV scheme : Ajjy = —1 Ffjy = zfnv = 0. (3-5) 
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In the HV scheme, vertex diagrams are treated differently in Ref. [|T 2 j and Ref. ||T5||. As 
a consequence, as already mentioned above, in the HV scheme one has to multiply the 
coefficients C±(iu) defined in refs. fl 0 |,[l 2 [] by the factor 1 — 4 a s (/x)/( 37 r), while this factor 
is already included in the definition of Ref. ]T5|]. We checked that our final result is 


identical in the NDR and HV schemes up to terms beyond NLO accuracy, if we use the 
expressions for C±(ji) of Sect. |II| including the additional factor just mentioned in the 
HV scheme. 

The coefficient functions quoted in Sect. [H] refer to a Fierz version of the weak Hamil¬ 
tonian different from (2T) and Fierz transformations do not commute with renormaliza¬ 
tion in general. If we use the standard Fierz version rather than the singlet-octet form 
quoted in ( |2.1|) , this interchanges Y R and Z R in the results quoted in Appendix A.3. How¬ 
ever, since in both schemes we used one has Y R = Z R , either of the two Fierz versions 
can be used. 

The NLO calculation for n = 3 5'( li has already been done in Ref. |J in the HV scheme. 
We find that our result for the functions gj[ 3 S'{ 1) ] defined in ( [2.13 ) and given in the 
Appendix agrees with the result of Ref. || . Nevertheless our result for the contribution 
of this channel to the decay rate, given by /[ 3 Si ](??) +NLO terms, differs from the 
one given in Ref. 0, because the authors of Ref. || used the coefficient functions of 
Ref. |L| , but did not correct them (or alternatively, the low energy matrix elements) by 
the factor 1 — 4 a s (/i)/( 37 r). As explained above with the conventions of Ref. [ 12 [ this 


additional factor is necessary in the HV scheme to obtain a scheme-independent result. 


B. IR regularization and NRQCD factorization 


The real and virtual corrections individually have double-logarithmic IR divergences, 
which we regulate by a gluon mass. However, after adding all contributions to the 
partonic process b —> cc+A", the IR divergences do not cancel completely. The remaining 
IR divergences are associated only with emission from the c and c quark. This is a 
necessary (but not sufficient) requirement for their factorization into NRQCD matrix 
elements as discussed in detail in Ref. [|l|. 

In addition to these IR divergences related to soft gluon emission, the last diagram in 
Fig. |I| exhibits the well-known Coulomb divergence, when the relative momentum of the 
c and c is set to zero. We regularize this divergence by keeping the relative momentum 
finite in the integrals, which would otherwise give rise to the Coulomb singularity. 

In order to extract the short-distance parts T[m] (see (|2.13|) ) of the partonic decay, 
we write 


r(6 -> cc[n] + X) = £r[m] ( O pert [m ]), (3.6) 

m 

where {O pert \m\) denotes the NRQCD matrix element for a perturbative cc pair in the 
state n. At NLO one has to calculate the left-hand side and ((9 pert [m]) to NLO. 






















FIG. 3. Perturbative corrections to the NRQCD operator matrix elements. A shaded circle 
denotes an insertion that specifies the angular momentum and colour of the cc pair. The 
vertical line implies that the diagram is ‘cut’. Symmetric diagrams are not shown. 


The diagrams that contribute the a s correction to ((9 pert [m]) are shown in Fig. [3]. 
For the first diagram (together with its complex conjugate) we obtain 


(O pert ' 1 [n]) A = (O pert ’°[n]) • A[n] a s 


2n 2 


(3.7) 


where A[n\ = Cp = 4/3, if n is a colour singlet state, A[n] = —1/(2 N c ) = —1/6 if n is a 
colour octet state and v is the relative velocity of the two quarks. (The superscript ‘O’ 
refers to a matrix element at tree level, ‘1’ denotes a 1-loop contribution.) This renders 
the short-distance coefficients free of the Coulomb singularity. 

The other two diagrams (called collectively ‘B’) together with their symmetry part¬ 
ners are UV and IR divergent. We define the NRQCD matrix elements in the MS scheme 
and denote their renormalization scale by fi. The IR divergence is regulated with a gluon 
mass to be consistent with the IR regulator used for the evaluation of the partonic pro¬ 
cess on the left-hand side of (|3.6| ) . The result is (compare with the Appendix of Ref. [[I] 
and with Ref. pAl, where other IR regulators are used): 


(or u ( s si)>B = ^ (in k+ T ( -48) (or-°( 3 a))K, 


(CT' I ( 3 Si))B = ? i (lnA + l 


A 2 


47T 




■f) (orwiK 


+ 


-20) <or*’YPo)>/ml 

(or u (‘so))B = 

(or u (‘so))B = 

+ (-f) (orh'n))/™? 


a s 

/, A 2 1 



An 

V R 2 3 

a s 

/ A 2 1 


In — — 

An 

V R 2 3 


16 

T 


32 

27 


(or tfi CPi))/m 


(3.8) 
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(A denotes the gluon mass.) Note that if one breaks up the 3 P term into terms with 
different J, one should replace 

<Or?'°(3po)} ^ 1/9 ((0 p, rt ,o ( 3 Po)) + (O Pjt.0 ( 3 Pi)) + ( O P|t.0 (3 p 2 ))) (3. 9) 

Using these results and solving for T[n] we hnd the IR hnite short-distance coefficients 
for each n collected in Appendix A.3. 

C. Difficulties with the colour singlet channels 

The LO contributions to the colour singlet channels 3 S[^ and A P[ 1 ^ are proportional 
to the small and strongly scale dependent coefficient (7m (//). One would therefore expect 
the NLO contribution to be particularly important for these channels. However, the 
strict NLO calculation leads to a negative, and therefore meaningless decay rate into 
these channels and to the conclusion that a reliable result can only be obtained at next- 
to-next-to-leading order. This problem was already identified and discussed in Ref. 0. 
(For the remainder of this section it is assumed that the reader has consulted Ref. [0 
for more details.) 

Consider the three next-to-leading order terms g t in ( |2.13| ). Despite its large coeffi¬ 
cient (7j|] the ( 73 -term, which comes only from a real correction, turns out to be numeri¬ 
cally very small (see the tables in the following section). Both g x and g 2 (at /i = are 
large and negative, and g 2 in particular, which comes with the larger coefficient 2(7[i](7[8], 
drives the decay rate negative. 

The authors of Ref. [[J suggested treating the decay process in a simultaneous expan¬ 
sion in a s and Upj/Cjg]. This implies that one should add to the term of order a s C^C\^ 
all terms of order a 2 C'L, because they also count as NLO in this rearranged expansion. 
On the other hand, the term a s Cfn (which involves g x ) should be neglected as being of 
higher order. The authors of Ref. || did not actually calculate all terms of order a^(7r|i, 
but estimated them by adding 

r o<°>]>(^) 2 U>i7jf P- 10 ) 

to (|2.13|) . This estimate can be motivated as follows: the virtual contribution to g 2 is 
given by the first four diagrams in Fig. |I] times the (complex conjugated) tree amplitude. 
All two-particle [cc]q cuts to the a^(7j|j term are given by the square of the first four 
diagrams in Fig. [T|. Hence, ignoring the real contribution, one may argue that g 2 [n] 2 /f[n\ 
is close (but not equal) to the two-particle contributions to the a 2 C 2 8 ^ term. 

In Ref. ]T6| the square of the 1-loop amplitude with a [cc]g final state is computed 
exactly and argued to provide a better estimate than the original one of Ref. || , because 
one leaves out only real contributions to the coefficient of a 2 C'L, which are argued to 

be phase-space suppressed. However, we hnd that for the :i P x 1 channel the virtual 
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contributions alone are IR divergent. Therefore the real correction that cancels this 
divergence cannot be argued to be small. In our opinion this also calls into question the 
assumption that the real contributions are numerically small for the S'-wave channels. 
For this reason we choose to follow with a minor modification the procedure of Ref. |J, 
which adds an IR finite term by construction, since g 2 is IR finite. The minor modification 
is the following: the third and fourth diagrams in Fig. [T] have imaginary parts, which 
contribute to the real part of the square of the amplitude (and hence to the coefficient of 
a 2 s C‘L ). The remnants of these imaginary parts after multiplying the one-loop amplitude 
by the complex conjugate of the tree amplitude can easily be restored from the ln(l — rj)- 
term (and Li 2 (? 7 )-terni in the case of ^P\) in the results presented in Appendix A.3. 
If we call the restored imaginary part Img 2 , then we use ( R10|) , with g% replaced by 
gl + (Img 2 ) 2 . 

We wish to emphasize two points: first, the discussed modifications of the colour 
singlet channels are certainly ad hoc and should be regarded with great caution. Second, 
the effect on the decay rate into a particular quarkonium state is not severely affected 
by this uncertainty, because it is dominated by colour octet contributions, whose short- 
distance coefficients can be computed reliably at NLO as we shall see. 

In order to gain a numerical understanding of the importance of the various terms 
involved in the colour singlet channels, we consider in the following three computational 
schemes for the decay rate: (a) the (strict) NLO calculation; (b) the NLO calculation 
with the term (|3.10|) added, but without the < 71 -term (‘improved’); (c) the same as (b), 
but with gi included (‘total’). For the S -wave colour singlet channels, (a) and (c) yield 
a negative rate. They are therefore meaningless. Option (b) yields a positive result of 
a magnitude similar to the result of Refs. [||[TB|. It may be considered as an order-of- 
magnitude estimate for the colour singlet contribution, but it may well be uncertain by 
100%. For the 3 P 1 - 1 ^ channel all three options give negative partial rates. However, since 
this channel mixes with and only the sum of the two is physical, a negative partial 
rate is not unphysical by itself. 


IV. RESULTS AND DISCUSSION 


In this section we present our results for the branching fractions of B decay into char- 
monium and moments of the quarkonium energy distributions in numerical form. The 
analytic expressions that enter ( |2.13| ) and (|2.15|) are collected in the appendices for ref¬ 
erence. 


A. Branching ratios for B decay into charmonium 
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1. General discussion of NLO corrections 


We normalize our calculation to the theoretical expression for the inclusive semileptonic 
decay rate 


r ,h _ Gl\V^mi 

1 SL “ 


192vr 3 


(l - 8z 2 + 8z 6 - z 8 - 24 z 4 Inz) m (z), 


where z = m c /m b and 


ffiO) = 1 


2 a s (m b ) {3 
3n 


+ ( ~~~T + 71-2 ) (!--) 


(4.1) 


(4.2) 


represents an excellent approximation |T7| for the 1-loop QCD correction factor. (The 
complete analytic result can be found in Ref. fl8|.) For any particular quarkonium state 
H , we obtain the branching fraction in the forndj 


Br ( b^h + x )= mj 2 (o"N) [ql 8] /NM (i + M«l) 

n 

[pfi]9i[ n }( 7 l) + 2C[i]C'[ 8 ]5f 2 M(^) + C 2 8 \g 3 [n)(ri)^ 


+ 


(4.3) 


The overall factor is given by 


u = Br“ p -5- = 3.0 • 10” 2 GeV” 3 , (4.4) 

^sl 

where we used BrgL P = 10.4% and To as given by ( |2.14|) . The charm and bottom pole 
masses are taken to be 1.5 GeV and 4.8 GeV, respectively. This yields 7] = 0.39, which 
we use unless otherwise mentioned. The sensitivity of the charmonium production cross 
sections to the quark mass values will be discussed below. 

We first examine the impact of the next-to-leading order correction and the depen¬ 
dence on the factorization scale /i for each intermediate cc state separately. We neglect 
the penguin contribution for this purpose. In Table | we show the branching fractions 
excluding the dimensionless normalization factor Af(0 H [n]) for three values of /j at LO 
and at NLO. To evaluate the LO expression we also use the Wilson coefficients at LO 
and 1-loop running of the strong coupling with A q^ d such that a s (M z ) = 0.119 both in 
LO and NLO. This is a large effect for the colour singlet channel, since Corrib) = 0.55, 
C[s](m&) = 2.14 (in the NDR scheme) but C^{m b ) = 0.41, C^fm b ) = 2.19. 


'When n is a P -wave state, (0 H [n\) should be understood as (0 H [n])/mf in the following 
formula, so that all matrix elements have mass dimension 3. Furthermore, in the case of 
n = which refers to the cc state 3 p( 8 ) summed over J = 0,1,2, the NRQCD matrix 

element (0 H [n\) is chosen to be (Og ( 3 Po))/m 2 . 
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n 

LO 

NLO 

/V GeV 

2.5 

5 

10 

2.5 

5 

10 

K’ 

0.0453 

0.201 

0.407 

-0.219 

-0.426 

-0.554 


0.0269 

0.119 

0.242 

-0.119 

-0.250 

-0.334 

3 C„ (1) 

0 

0 

0 

-0.660 

-0.481 

-0.377 

3 ipf 1} 

0.0537 

0.238 

0.484 

-0.654 

-0.738 

-0.794 

3 p(l) 

0 

0 

0 

-0.534 

-0.389 

-0.305 

lpW 

0 

0 

0 

-1.02 

-0.741 

-0.58 

i s (s) 

8.72 

8.01 

7.51 

12.6 

11.1 

10.2 

3^(8) 

5.18 

4.75 

4.46 

7.70 

6.80 

6.18 

3 Pj 8) 

31.1 

28.5 

26.8 

38.3 

34.5 

31.7 

1 P{ 8) 

0 

0 

0 

-1.95 

-1.53 

-1.30 


TABLE I. Comparison of LO and NLO for the decay rate into a cc pair in state n and 
the dependence on the factorization scale p. The dimensionless overall factor A f (0 H [n\) is not 
included. Quark masses: mj, = 4.8 GeV, m c = 1.5 GeV. Penguin contribution not included. 


We now observe that the colour singlet contributions are, as expected, enormously 
scale-dependent at LO. The NRQCD matrix element (0{^( 3 Si)) is related to the radial 
wavefunction at the origin by (O'^^( 3 Si)) = 9|i?(0)| 2 /(27r) up to corrections of order u 4 . 
Using (Oi^( 3 Si)) = 1.16 GeV 3 |I9| , we obtain 


Br (B —> J/ip + X ) = (0.09 — 0.84)% (colour singlet, LO) 


(4.5) 


to be compared with the measured branching fraction (0.80 ± 0.08)% BE The LO 
prediction is uncertain by a factor of about 10 for all colour singlet channels, as can be 
seen from Table |. As also seen from this table, the scale uncertainty is reduced to a 
factor 2-3 at NLO. However, the NLO correction term renders the partial decay rates 
negative, as already mentioned in Section |H1 C| . 

The situation can be somewhat improved by adding the estimate ( |3.10| ) for the order 
cc 2 NNLO term with the large coefficient CL, while treating the a s C r 2 ! term as formally 


iNote that we denote by Br (B —> J/'ip + X) the direct production of J/ip, excluding radiative 
decays into J/ip from higher-mass charmonium states. The same convention applies to all other 
charmonium states H. 
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n 

/i 

W 

3 S W 

spa) 


2.5 

0.0453 

0.0269 

0.054 

LO 

5 

0.201 

0.119 

0.238 


10 

0.407 

0.242 

0.484 


2.5 

-0.219 

-0.119 

-0.654 

NLO 

5 

-0.426 

-0.250 

-0.738 


10 

-0.554 

-0.334 

-0.794 


2.5 

0.0951 

0.034 

-0.392 

‘Impr’ 

5 

0.058 

0.0259 

-0.280 


10 

0.0706 

0.0402 

-0.176 


2.5 

0.0283 

-0.0138 

-0.455 

‘Tot’ 

5 

-0.0984 

-0.0862 

-0.427 


10 

-0.159 

-0.124 

-0.391 


TABLE II. The same quantity as in Table I for various treatments of the colour singlet 
channel (for which the LO short distance coefficient, i.e. f[n], is non-vanishing). Penguin 
contribution not included. 

of higher order in a double expansion in a s and Cpj/C'jg] @. The addition of ( p.lOj ) 
also reduces the factorization scale dependence further, because it contains exactly the 
double logarithmic correction ln 2 (m^//x 2 ), which is required to cancel the large 

scale dependence of at leading order in a s . In Table |TT] we display the result for 
the partial decay rates into the colour singlet channel, which is obtained in this way 
(denoted ‘Impr’ in the table) and for comparison again the LO and NLO result. The 
improvement can and should be done only for those colour singlet channels that have 
non-vanishing LO contributions. The last three rows of Table [H] show the results that 
are obtained if we add back the g\ term in ( |4.3D to the improved treatment. The a s C'm.9i 
term is sizeable and negative and therefore re-introduces a large scale dependence. The 
same improvement that is applied to the LO C^f term is necessary for the gi term, 
which would require going to order ot^Cf s y One may argue that unless this is done, it is 
preferable to leave the g\ term out entirely. Therefore we shall use the ‘improved’ version 
(‘Impr’ in Table p) as our default option later. While the result is certainly not accurate, 
we believe that this is the best we can do to the colour singlet channel without making 
arbitrary modifications. We note that for J/^ this gives a colour singlet contribution to 
the branching fraction, which is close to the lower limit in (|4.5|) and also compatible with 
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the estimates of Refs. [:9j.[L6|. ft seems safe to conclude that colour singlet production 
alone is not sufficient to explain the measured branching fraction. 

The partial rates in the four relevant colour octet channels are shown in the lower 
part of Table ||. In this case we fold that the perturbative expansion is very well behaved. 
The NLO short-distance coefficients are larger by 20%±50% than the LO coefficients and 
the scale dependence is very moderate. The scale dependence is not reduced from LO to 
NLO. This is due to the fact that the LO coefficients depend only on the scale-insensitive 
C[g], while there are sizeable coefficients of the highly scale-dependent combinations 
C'jijC'jg] and Cm at NLO. The numerical enhancement of the short-distance coefficients 
in the colour octet channels, which is evident from Table is sufficient to account for 
the measured J/0 branching fraction, as already noted in Refs. []7]||. The positive NLO 
correction reinforces this trend. Other J/0 production processes suggest that the long¬ 
distance parameters in the colour octet channels are of the order a few times ICO 2 GeV 3 . 
(This will be made more precise soon.) This leads to typical branching fractions of order 
0.5%. 

For our standard value of the quark masses (= 4.8 GeV, m c = 1.5 GeV), the 
numerical values for the functions that enter (|4.8| ) are given in Tables pTIj - |Vfff| for all of 
the six (known) charmonium states below the DD threshold. We keep the dependence on 
the factorization scale of the weak Hamiltonian /i, but put the NRQCD factorization scale 
equal to 2 m c . At the scale = 4.8 GeV we also have a s (rrib ) = 0.22 and Corrib) = 0.55, 
C[ 8 ](m 6 ) = 2.14. The tables also include the penguin correction factor 5 P . 

We fold that the dependence on the quark masses changes little when going from LO 
to NLO (for those n for which the LO term is non-zero). The quark mass dependence is 
reasonably well estimated by that of the ratio 


(1 — 4m 2 /m 2 ) 2 
mlm c fi(m c /m b ) ’ 


(4.6) 


where fi(z) is the tree-level phase space factor in round brackets in (3~I) and the NRQCD 
matrix elements are assumed fixed. If we vary m c by ±100 MeV and mb by ±200 MeV 
around our ‘standard’ values and add the separate variations in square, we find a variation 
of r of about 15% relative to the standard value. Taking into account the approxima¬ 
tiveness of this estimate, we assign an overall normalization uncertainty of 20% due to 
quark masses. 


2. B —> 0(nS’) ± X 

We now turn to a more specific discussion of B decay into the spin-triplet S'-wave states 
J/0 and 0', denoted collectively by 0( nS ) or just 0. For quick reference we give the 
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n 

f 

5p 

9i 

92 

93 

3 sf } 

0.661 

-0.004 

-20.5 

-8.46 + 2.65 ln(mg//z 2 ) 

0.162 

3^(8) 

0.992 

-0.09 

0.486 

-6.35 + 1.98 \w{ml/n 2 ) 

32.0- 3.97 \n{ml/n 2 ) 

lg(8) 

1.67 

0.009 

0.556 

—11.2 + 3.34 In {ml/ n 2 ) 

50.2 - 6.68 ln(m 2 b /ii 2 ) 

3 Pf 

5.95 

0.009 

-129 

-36.6 +11.9 In (ml/n 2 ) 

121-23.81n [mlhi 2 ) 


TABLE III. Numerical values for the LO and NLO functions (NDR scheme) that enter the 
branching ratio Br (B —> J/ip + X) according to (A3) for mi, = 4.8GeV, m c = 1.5GeV. The es¬ 
timate for the penguin correction 5p is obtained with the parameters detailed in Appendix A.2. 
The table applies without modification to ij/. 


n 

f 

5 P 

9i 

92 

93 

1 cl 1 ) 

1.11 

0.07 

-28.6 

—15.0 + 4.46 In {ml//i 2 ) 

0.185 

00 

1.67 

0.009 

0.556 

-11.2 + 3.34 In (m 2 b /n 2 ) 

50.2 - 6.68 ln(mg//x 2 ) 

CO 

'oo 

0.992 

-0.09 

0.486 

-6.35 + 1.98 ln(m 2 b /ii 2 ) 

32.0 - 3.97 ln(m 2 b /ii 2 ) 

lpM 

0 

0 

-28.2 

0 

-17.4 


TABLE IV. Numerical values for the LO and NLO functions (NDR scheme) that enter the 
branching ratio Br (B —> r] c + X) according to (4.3) for = 4.8 GeV, m c = 1.5 GeV. The esti¬ 
mate for the penguin correction dp is obtained with the parameters detailed in Appendix A.2. 


n 

f 

dp 

9\ 

92 

93 

3 C„ (1) 

0 

0 

0 

0 

-6.11 

3 

0.992 

-0.09 

0.486 

-6.35 + 1.981n (jnUfi 2 ) 

32.0 - 3.97 In (m 2 b //i 2 ) 


TABLE V. Numerical values for the LO and NLO functions (NDR scheme) that enter the 
branching ratio Br (B —► XcO + X) according to (4.3) for m& = 4.8 GeV, m c = 1.5 GeV. The es¬ 
timate for the penguin correction 5p is obtained with the parameters detailed in Appendix A.2. 
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n 

f 

5p 

9i 

92 

93 

3 P 1 (1) 

1.32 

0.07 

-26.8 

-16.3 + 5.291n (m 2 b /ii 2 ) 

-4.04 

3 s[ 8) 

0.992 

-0.09 

0.486 

-6.35 + 1.98 \n(ml/n 2 ) 

32.0- 3.97 \n(m 2 b /n 2 ) 


TABLE VI. Numerical values for the LO and NLO functions (NDR scheme) that enter 
the branching ratio Br (B —> Xci + V) according to (L3) for mj = 4.8 GeV, m c = 1.5 GeV. 
The estimate for the penguin correction 5p is obtained with the parameters detailed in Ap¬ 
pendix A. 2. 


n 

f 

5p 

9i 

92 

93 

CO 

fcO 

0 

0 

0 

0 

-4.94 

3 s[ 8 ^ 

0.992 

-0.09 

0.486 

-6.35 + 1.981n (jn\hi 2 ) 

32.0- 3.97 \n(ml/n 2 ) 


TABLE VII. Numerical values for the LO and NLO functions (NDR scheme) that enter 
the branching ratio Br (B —> \ci + X) according to (L3) for mj = 4.8 GeV, m c = 1.5 GeV. 
The estimate for the penguin correction 5p is obtained with the parameters detailed in Ap¬ 
pendix A. 2. 


n 

/ 

5p 

9i 

92 

93 

i p W 

0 

0 

0 

0 

-9.41 

1 c( 8 ) 
‘-’O 

1.67 

0.009 

0.556 

-11.2 + 3.34 In (ml/fi 2 ) 

50.2 - 6.68 In {ml/fi 2 ) 


TABLE VIII. Numerical values for the LO and NLO functions (NDR scheme) that enter 
the branching ratio Br (B —> h c + X ) according to (4T) for mb = 4.8 GeV, m c = 1.5 GeV. 
The estimate for the penguin correction dp is obtained with the parameters detailed in Ap¬ 
pendix A. 2. 
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branching ratio in completely numerical form for p = mj, = 4.8 GeV and m c = 1.5 GeV:| 


'-0.7411 

Br ( B -> iP(nS) + X) = <j 0.0754 

-0.254 


> 10~ 2 (Of ( 3 ^i)) + 0.195((9g ( 3 Si)) 


(4.7) 


0.342 


(of (‘So)) + yy (£>%(%)) 


3.10 


^/3r 


m; 


The penguin correction is included. For the coefficient of the colour singlet operator we 
display the result obtained according to the procedures ‘NLO’, ‘Improved’ and ‘Total’ 
in Table [IT]. As discussed earlier, we use the second entry (‘Improved’) in the following. 
The colour singlet matrix element is computed from the wave functions at the origin 


obtained with the Buchmiiller-Tye potential as given in Ref. |19[1: 


(Of( 3 S i)> = 


9|fl(0)| ! 

27T 


1.16 GeV 3 
0.76 GeV 3 


(W) 


(4.8) 


The colour octet matrix element (Og( 3 Si)) is rather well determined by direct 7p pro 
duction at large transverse momentum in pp collisions [[20] . We use the values |22 


(oh 3 Si)) 


1.06-10- 2 GeV 3 (J/ip) 

0.44 ■ 10- 2 GeV 3 {p/) 


(4.9) 


There is an uncertainty of a factor 2 in each direction of the central value associated 
with these numbers. With the number quoted the :i S[ S> channel contributes 0.21% to the 
J/p) branching fraction and 0.09% to the p;' branching fraction. The other two colour 
octet matrix elements are not yet well determined. In Fig. (|we show the J/pj branching 
fraction as a function of the renormalization scale p for various values of M^ A ( 1 5'{ 8 ' 1 , 3 Pj 8 ^), 
where 

Mp'sP, 3 Pj 81 ) = (0*(‘S „)) + A (Oi( 3 P 0 )). (4.10) 

m c 

With the other parameters fixed, the branching ratios Br (B —> J/p> + X ) = (0.80 ± 
0.08)% and Br (B —> ip' + X) = (0.34 ± 0.05)% measured by CLEO are reproduced by 


Mt ii 


‘Sf 


(8) 3p(8)> 


J 


1.5-lO” 2 GeV 3 (J/ip) 

0.6-10- 2 GeV 3 (ip'). 


(4.11) 


^In this equation and those of the following ones that are 
given in units of GeV~ 3 . 


similar in form, all numbers are 
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3456789 10 
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FIG. 4. The J/ijj branching fraction as a function of renormalization scale p for various 
values of Mf 1 ( 1 S[ S \ 3 P < f' > ) in 10“ 2 GeV 3 . The horizontal band shows the CLEO measurement 
[2] and the dashed curve the colour singlet contribution alone. 


If we allow the colour singlet contribution to vary between zero and twice the value 
assumed above, include the above variation of (Og( 3 Si)) as well as the experimental 
error, and add all variations linearly, we obtain the allowed range 


, 3 p( sh 


j 


(0.4 — 2.3) ■ 10 ~ 2 GeV 3 {J/ip) 

( 0.0 - 1 . 0 ) • 10 - 2 GeV 3 (?//)• 


(4.12) 


It is interesting to compare the central values ([4.11 ) and the upper limits with other 
determinations of the parameter ( 1 S[ 8 \ 3 Pj 8 ' > ). The central values are about a factor 
3 smaller than the central values obtained for M^’ 5 ( 1 S[ 8 \ 3 Pj S ' > ) from ^ production at 

m 


the Tevatron at moderate transverse momentum 


As emphasized in Ref. 


the Tevatron extraction is very sensitive to various effects that affect the transverse 
momentum distribution. Indeed, Refs. fl23l,f24| quote smaller values compatible with, or 


smaller than the central values above. The total production cross section in fixed target 
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collisions probes Mf ( 1 S[ 8 \ 3 Pj 8 ' > ) (assuming the validity of NRQCD factorization, which 
may be controversial). Given that a different combination of matrix elements enters, 
the values obtained in Ref. p5| are certainly consistent with the above central value. 


In view of the uncertainties involved in charmonium production in hadron collisions, we 
believe that the above upper limit on , 3 Pj 8 ^) is the most stringent one existing 

at present. We note that small values of (Of ( 1 S' 0 )) and (Og( 3 P 0 )) seem to be preferred 
by the non-observation of a significant colour octet contribution in the energy spectrum 
of inelastic J/-0 photoproduction |[2Ci|.[T|.|27| (see, however, the discussion in Ref. |l3| ). 


We conclude that the measured J/ip and ip’ branching fractions can be accounted for 
with values of the NRQCD long-distance parameters consistent with previously available 
values. 


3. B —»■ ? 7 C + X 


Presently, only an experimental upper bound Br (B —> 7/ c + X ) < 0.9% @] exists on r/ c 
production. For the same choice of input parameters as above, we have 


Br (B —> rj c + A") 



> IQ’ 2 {O^CSq)) + 0.342(OJ c ( 1 5o)) 


(4.13) 


+ 0.195 


(ores,)} 


0.240 

ml 


(Or('P,)) 


The LO term is enhanced by about 10% because of the penguin correction. 

There is at present no information on the r] c colour octet matrix elements from other 
7 / c production processes. The colour octet matrix elements are non-zero because soft 
gluon emission connects the colour octet cc state to the physical charmonium state. The 
soft gluon emission amplitude can be multipole expanded, supposing that the character¬ 
istic momentum of the emitted gluons is of order m c v 2 , smaller than the characteristic 
momentum m c v of the charm quarks in the charmonium rest frame. Up to corrections of 
order v 2 , spin symmetry imposes relations between the r] c and J/^ matrix elements. In 
addition to the familiar spin symmetry relation {0'^ u ’( 3 Si)) = 3 (O’^^So)) for the colour 
singlet wave function, we find 

(ores,))= i (o^es,)), 

(ores,)) = (o^eso)), ( 4 + 4 ) 

(orep,)} = 3 (o j e*ep 0 )). 


Note that these relations are consistent with the renormalization group equations for 
the matrix elements that follow from (3J3). Since we do not know ((9g // ^( 1 S'o)) and 
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(£>J/V’( 3 p o )) se p ara t e ly 5 we assume that between one half and all of Myp ((S'®, 3 P®) 
due to (O^fSo)). With the central value from ( |4.11|) this leads to the estimate 


is 


Br (P —> i~i c + X ) ps (0.3 — 0.5)%. 


(4.15) 


We emphasize that this estimate is crude and depends sensitively on the validity of the 
relations ( |4.14| ) . This estimate is below the J/ip branching fraction, but with the increase 
in statistics since the previous analysis [3|, a branching fraction in the above range may 
perhaps be reached with the CLEO detector. 


4- B —* Xcj + X 


Colour octet effects in charmonium production were in fact considered for the first time 
for x c production in P decay |5|. The authors showed that the observed y c2 signal can 
be explained by the production of a ccpP®] state followed by a soft dipole transition. 
(Recall that at LO in the colour singlet model, XcO and Xc 2 are not produced.) The 
LO production through a ceps'®] pair corresponds to the IR sensitive contribution at 
order a s in the ordinary colour singlet channel. Our NLO calculation adds to those the 
‘hard’ contributions at order a s in the colour singlet and the colour octet channels. With 
rrii, = 4.8 GeV, m c = 1.5 GeV and // = mb as usual we obtain 

Br (P - Xc0 + X) = ^°^ 2 148 (Of°( 3 P 0 )) + 0.195 (0^%%)), (4.16) 

in -2 

—y (O^CPi)) + 0.195 (0 8 Xc 1 ( 3 S!)), (4.17) 

m 2 

Br (B ^ Xc2 + X) = ^ 0 -° 2 120 (Of 2 ( 3 P 2 )) + 0.195 (0^ 2 (%)) (4.18) 

to be compared with the measurements [| 2 ] 


Br (P —> Xd + X) = l -0.783 

- 1.21 


Br (P Xci + X) = (0.37 ± 0.07)%, (4.19) 

Br (P -»■ Xc 2 + X) = (0.25 ± 0.10)%. (4.20) 

Owing to the spin symmetry relations, valid up to higher order corrections in v 2 , 

(Oi cJ ( 3 Pj)) = (2 J + 1 ) (Of°( 3 P 0 )), (4.21) 

(or- / ( 3 s 1 )) = (2 J + 1 ) (or 0 ^)), (4.22) 


only two of the six above parameters are independent. The colour singlet contribution is 
always negative. However, since the two matrix elements involved mix under renormal¬ 
ization each short-distance coefficient depends on an arbitrary convention to separate 


21 











the two contributions. (We used the MS scheme.) Hence a negative partial rate is not 
unphysical. 

Due to the near-proportionality of ( |4.17|) , (|4.18|) , we find, however, that it is very 
difficult to reproduce the experimental result with a reasonable choice of matrix elements. 
With most reasonable choices one obtains a \ci production cross section larger than the 
cross section for Xci ■ If we take the colour singlet matrix element from the Buchmiillcr- 
Tye potential model (JHJ, 


(Of c0 ( 3 P 0 ))/m 2 c = 4.8 • 1(T 2 GeV 3 , 

and adjust^ 

(C>f 0 ( 3 Si)) = (4.5 - 6.5) • 1(T 3 GeV 3 
to reproduce the measured Xc 2 branching fraction, we obtain 

Br (B -> Xd + X ) = (0.15 - 0.27)%, 


(4.23) 


(4.24) 


(4.25) 


which is below the measurement. We conclude that the expansion in a s is not well 
behaved enough for P -wave charmonium production in B decay, at least to next-to- 
leading order, to arrive at quantitative relations. The problem is caused to a large 
extent by the fact that the leading order colour singlet contribution to Xci production, 
which would have been expected to enhance Xci production relative to Xci production, 
is turned negative (or very small) at NLO. 

In addition, the colour singlet contributions are also negative for Xco and Xc 2 , which 
have no leading order contribution. This requires large cancellations between a negative 
colour singlet contribution and a positive colour octet contribution. These cancellations 
may be considered an artefact of the MS factorization scheme, which appears unnatural 
from this point of view. 

Because of this unsatisfactory situation, we fold it difficult to predict the Xco branch¬ 
ing fraction better than the naive expectation of one fifth of the Xc 2 branching fraction. 
Note that we consider the prediction for the Xci state less reliable than for the % c2 state, 


so that the range for the colour octet matrix element obtained in (|4.24|) may not be 
totally unreasonable. 


5. B —> h c + X 

Finally, we consider the production of the l P\ state h c , which will probably be difficult 
to measure. With the usual choice of parameters: 


**This value can be compared with (Og c0 ( 3 S i)) = 3.2 • 10 3 
Xcj production at the Tevatron collider. 


GeV 3 obtained in Ref. [21] from 
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Br (B^h c + X ) = 0, ° 2 228 (O^CPi)) + 0-342 (0^(X)). (4.26) 

m 2 c 

We use the spin symmetry relation (0\ C { } P\)) = 3 ((9* c 0 ( 3 Po)) for the derivative of the 
wave function at the origin (squared) and the spin symmetry relation 

(O h 8 c( 1 S 0 )) = 3(Or°( 3 S 1 )) (4.27) 

for the colour octet matrix element. With the value of the colour singlet matrix element 
as quoted above and the colour octet matrix element in the range ([4.24]) we obtain the 
estimate 


Br (B -> h c + X) « (0.13 - 0.34)%. 


(4.28) 


The given range hinges on the estimate (|4.27|) , which follows from the statistical factor 
3 for the El (% — 1 ► l P\ transition and another factor 3 for the final state with J = 1, 
provided we assume a spin-independent overlap integral in terms of colour singlet and 
colour octet wave functions. Then taking into account the fact that the definition of 
(Og^^Si)) does not include an average over the three initial polarizations, we obtain 
the relative factor 9/3 = 3. Note again that (|4.27| ) is consistent with (|3.8|) . 


B. Quarkonium momentum distributions 


We briefly address quarkonium momentum distributions. We restrict ourselves to the 
S-wave states J /-0 and //, Erst because data exist at present only for these states [Q and 
second, because the theoretical prediction appears to be most reliable for the S'-wave 
states. 

The energy distributions of the cc -pair in the partonic process b —> cc[n] + q + g are 
collected in Appendix B. The partonic energy/momentum distributions are distributions 
in the mathematical sense and cannot be used to predict the physical momentum spec¬ 
trum. Two effects lead to a smearing of the partonic energy/momentum distribution: 

(a) The b quark has a residual motion in the B meson and does not decay at rest. 
This leads to an energy smearing of order A, where A is the QCD scale. For B decay 
into charmonium this ‘Fermi motion’ effect has already been modelled in Ref. [[28[ 


and 


more recently in Ref. [29 


(b) The second effect is related to the fact that charmonium production through 
colour octet cc[n) states requires the emission of soft gluons with energy of order 
m c v 2 ~ A. The NRQCD matrix elements measure the probability for the transition 
n —> charmonium, but do not take into account the kinematic effect of the soft gluon 
emission, which becomes important near the upper endpoint of the charmonium momen¬ 
tum distribution 0- Because the maximal J/'ij) momentum in B decay is only about 
2 GeV, this effect is expected to be at least as important as the b quark motion, especially 
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Moment 

(Of( 3 so> 

(O s *( 3 Si)) 

(O s *(B„)> 

( 0 *( 3 P„))/ml 

0 

7.5 • 1 (T 4 

2.0 • 10” 1 

3 . 4 - 10” 1 

1.06 

LI 

5.1 • 10~ 4 

1 . 9 - 10- 1 

3 . 3 - 10” 1 

9.8 • 10” 1 

L 2 

3.3 • 1 (T 4 

1 . 8 - 10- 1 

3.2 • lO ” 1 

9.2 • 10” 1 

L 3 

1 . 7 - 10 -4 

1 . 7 - IQ ” 1 

3.1 ■ 10” 1 

8 . 7 - lO ” 1 

SI 

2.4 • 1 (T 4 

00 

o 

1 

CO 

1 . 4 - 10“ 2 

CD 

o 

l 

to 

S 2 

5.4 ■ 10“ 5 

1.6 ■ 10“ 3 

2.3 ■ 10 -3 

1.6 • 10“ 2 

S 3 

2.4 • 10“ 5 

5 . 8 - 10“ 4 

7.5 • 10" 4 

6.1 • 10“ 3 


TABLE IX. Coefficients of NRQCD matrix elements in GeV~ 3 for the moments of the 
J/'ij) momentum distribution. ‘Ln’ denotes the moment with the weight function z n , ‘Sn’ the 
moment with the weight function (1 — z) n . 


in view of the large fraction of colour octet production. No attempt has been made so 
far to take this effect into account in a realistic model for the spectrum. Both effects 
are also related to the fact that the phase space boundaries in the partonic calculation 
depend on the quark masses rather than the B meson, J/ijj and K masses. 

In this paper neither of the two effects will be modelled. One expects them to be 
suppressed by powers of A QCD/ m b,c and v 2 provided a smooth average of the momentum 
distribution (like the integration to the total width) is taken. We define 


Br {$, W) 


i 

J dzW(z) 

o 


1 dr{B^n(> + X) 
r B dz 


1+7/ 
f dx 


n 




dr [n] 
dx 


W 


f y/x 2 — 4rj\ 

V l-*7 )' 


(4.29) 


where z = \pip\/\p^,max\ and is the total B decay rate. The variable x is the energy 
fraction, defined together with dT[n]/dx in (|2.15|) . (The functions g t [n] (rj, x) that enter 
there can be found in Appendix B.) Finally W(z) is a smooth weight on the momentum 
distribution. 

In Table |IX| we give the coefficients in front of the NRQCD matrix elements in ( [f.29| ) 
for the weight functions z n (Ln) and (1 — z) n (Sn) up to n = 3. For n = 0 we recover 
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the inclusive branching fraction ( |4.7| ).p| The hrst weight function increasingly weights 
the endpoint region as n increases. Therefore one cannot take n large without enhancing 


higher order terms in the velocity expansion [13j not taken into account here. The second 
weight function weights the small-momentum tail. The moments decrease rapidly with 
n, because, as expected, the spectrum favours a large momentum of the ip. The LO order 
and NLO virtual contributions do not contribute to the Sn moments. These moments 
are directly sensitive to hard gluon radiation. 

The momentum spectrum measured by CLEO |2| is given in the CLEO rest frame 
rather than the B rest frame. With the improved statistics that should be available 
now, it will be interesting to see whether one can obtain additional information on 
charmonium production by comparing averages of the momentum spectrum with the 
above predictions. For example, one may think of using the Sn moments to determine 
the parameters (Og^So)) and (Og( 3 P 0 )) separately. 


C. Comparison with exclusive two-body modes 

It is interesting to compare qualitative features of the theoretical result for inclusive 
charmonium production with the sum of the most important exclusive two-body decay 
channels containing charmonium. We discuss only J/ip because only limited experimen¬ 
tal information exists for the other charmonium states. 

There exist only two two-body modes of any importance. Their branching fractions 
are measured to be |30 


Br (B -»■ J/ip ± K) = 
Br (B -> J/ip ± K *) = 


' (0.099 ±0.010)% 
(0.089 ±0.012)%, 

' (0.147 ±0.027)% 
(0.135 ±0.018)%, 


(4.30) 

(4.31) 


where the upper line in both (|4.30| ) and (|4.31|) refers to B ± and the lower one to B°. 
The combined branching fraction of about 0.25% is far below the inclusive branching 
fraction of 0.80 ± 0.08%. The existence of a large fraction of three-body modes is also 
confirmed by the broad energy distribution measured by CLEO @ . 

There is no clear association of the colour singlet contribution with the inclusive 
branching fraction with the above two-body modes. A colour singlet cc pair can radiate 


it The precise implementation is done as follows: For the delta-function term in ( p.l| ) we 
use the ‘improved’ prescription for the colour singlet channel. For the second term on the 
right-hand side of (B4) we use the strict NLO approximation. This is reasonable, because 
the negative contributions that necessitate the improvement are mainly associated with virtual 
corrections. 
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a hard gluon and end up as a three-body mode. Likewise, the soft gluon that is necessarily 
emitted from a cc pair in a colour octet state can be re-absorbed by the light quarks 
and combine to the kaon. This process depends on the spectator light quark in the B 
meson and would therefore seem to violate the factorization hypothesis of the NRQCD 
approach, ffowever, factorization requires only that such spectator dependence cancels 
out to first approximation (in A qcd/^c^) in the average over all decay modes. 

Nevertheless it is natural to expect that a cc pair in a colour octet state finds itself 
more often hadronizing in a multi-body decay than a colour singlet cc pair and there¬ 
fore we consider the observation of a large fraction of multi-body decays as supporting 
evidence for the colour octet picture. The fact that the sum of the two two-body modes 
above is larger than the (poorly predicted) inclusive colour singlet branching fraction of 
approximately 0.1% suggests that the NLO calculation underestimates this contribution. 
(Another possibility is that some fraction of the large colour octet partial rate does in 
fact end up in two-body modes.) 

The large fraction of three body decays is significant in another respect. The NRQCD 
approach assumes ‘local parton-hadron duality’, which is often discussed in connection 
with the inclusive non-lcptonic decays of B mesons. A crucial consequence of ‘local 
parton-hadron duality’ is that the effect of colour reconnections to the spectator quarks 
are small (power-suppressed) even though colour reconnections must occur in every single 
decay. The same assumption underlies the NRQCD factorization approach to inclusive 
charmonium production. The assumption is usually justified by the presence of a large 
energy release into the final state and many decay modes to be averaged over. Since 
the energy release in decays into charmonium is not particularly large, the existence of 
a sufficient fraction of decays with additional pions suggests that the total decay rate 
provides enough averaging for an approximate cancellation of non-faetorizable effects.^ 


V. CONCLUSION 

We presented an analysis of inclusive B decay into the known charmonium states at next- 
to-leading order in the strong coupling and accounting for the most important colour 
singlet and colour octet production mechanisms. We find that radiative corrections make 
the colour singlet contributions negative, an effect already observed in Ref. JSJ for the 
colour singlet contribution to J/-0 production and explained by the suppression of colour 
singlet production at leading order due to the particular structure of the weak effective 


■f+We use the term ‘non-factorizable’ in the sense of NRQCD factorization. In the literature on 
two-body decays of B mesons ‘non-factorizable’ usually refers to any virtual gluon correction 
(below the scale mb) that connects the b or s quark to the core quark. These ‘non-factorizable’ 
contributions are included in the inclusive NRQCD calculation, see Fig. [l]. 
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Hamiltonian. The problem is particularly serious for \ci production. As a consequence, 
the appealing theoretical pattern for XcJ production at leading order 0 becomes ob¬ 
scured quantitatively at next-to-leading order, although the qualitative requirement of 
an additional colour octet component is not put into question. In general, we find it 
difficult to make a quantitative prediction for any of the four P -wave states. 

The situation is more satisfactory for S-wave production, and in particular for J /if 
production, for which we confirm the earlier conclusion that the colour singlet contri¬ 
bution is about a factor 5-10 below the observed production rate. The next-to-leading 
order corrections to the colour octet channels computed in this paper are positive and of 
order 20%-50%. Assuming that these channels make up for the missing contribution, we 
adjust a certain combination of NRQCD matrix elements to reproduce the experimental 
branching fraction. The value for the matrix element is somewhat smaller than, but 
within errors compatible with the magnitude suggested by other quarkonium production 
processes. Since the a s expansion appears to be well-behaved for J/^ production, we ob¬ 
tain the sharpest upper bound to date on a combination of (O^ (^o)) and (0$ ( 3 P 0 ))/ml, 
even when we assume that there are no other production mechanisms. In principle, with 
more accurate data appearing, weights on the J/ft momentum distribution can also be 
used to determine (Og ( 1 i S'o)) and (Of ( 3 Po))/ml individually, which has so far proven to 
be difficult, even combining information from other J/tf production processes. 
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APPENDIX A. SHORT-DISTANCE COEFFICIENTS FOR INTEGRATED 

DECAY WIDTHS 


We collect the expressions that enter the integrated partial rates ( |2.13| ) in this appendix. 
Recall r] = 4ml/ml- The scale fi denotes the factorization scale at which the coefficients 
(7(1,8] are evaluated. We distinguish from it the renormalization scale ft, of the NRQCD 
matrix element (0 H [n\). 
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A.l. Lowest order functions 


The lowest order function f[n](r)) vanishes for n = 3 Pq 11 , Z P^\ { P ['^ and 'Ppx The 


non-vanishing ones are: 


/[ 3 ^ 1} ]('7) = (1 - '?) 2 (1 + 2/ 7 ), 

(A.l) 

SP 

GO 

1 

to 

(A.2) 

/pPf’lh) = 2(l- I? ) 2 (l + 2^). 

(A.3) 

/PsPlM = | (i -Pi 1 + 2 n), 

(A.4) 


(A.5) 

/pPf ](») = 9(l->?) 2 (l + 2u). 

(A.6) 

A. 2. The penguin correction 


We consider only the contribution where a decay through a QCD penguin operator 
03_6 interferes with the decay amplitude through a current-current operator. Since the 
penguin operators have small coefficient functions, it is sufficient to evaluate the penguin 
contributions in lowest orders in a s . The double penguin contribution is negligible in 
size. Since f[n](r)) = 0 for n = 3 Pq 1 \ ■, { P \'^> l P^\ the penguin contribution also 

vanishes in leading order for these n. For the other intermediate states, we find 

ipfSP] = 2 3(03 + + + Ce « 0.004. 

C[l] 

(A.7) 

= M 3 ?' 11 ] = 2 3(Cs - C =) + C * - C ‘ K 0.07. 

C[l] 

(A.8) 

Vpsf’l = 4 A + Ge « 0.09, 

C[ 8 ] 

(A.9) 

dpfs^} = 5p[ 3 pj 8) ] = 4 ° 4 ~ ° 6 ~ 0.009. 

C[ 8 ] 

(A.10) 


Here we used unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) matrix to relate 
the CKM factor of the penguin contribution to the CKM factor of the current-current 
operator contribution (up to a negligible error Re (V))) Ws/(b^Vt s )) = —1). For the nu¬ 
merical estimate we have used the leading logarithmic approximation for the Wilson co¬ 
efficients at the scale m b = 4.8 GeV and with A qc D adjusted to reproduce a s (mz ) which 
gives A Q£ D = 93MeV (for 5 flavours). This implies C[i](mb) = 0.41, C[ 8 ](mb ) = 2.19, 
C 3 (m b ) = 0.010, C 4 (m b ) = -0.024, C 5 (m b ) = 0.007 and C 6 (m b ) = -0.028. 
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A.3. Next-to-leading order functions 


The next-to-leading order functions depend on the scheme-dependent constants Xr, Yr 
and Zr. Their values in the NDR and HV scheme are given in section [111 A. 


A.3.1. ;: Sj l:i 


SipSf’lW = 3 (-l)(l->)) 2 (l + 2r,) 


81+2(77) + 41 n(l — rj) In77 — 4 Zr + 


47 T 


2 -i 


— 4 77(1 + 77) (1 — 2 rj) In 77 — 2 (1 — rj ) 2 (5 + Arj) ln(l — rj) 
-(I-77) (I-277) (3 + 577 )). 


(A.ll) 


92 [ 3 S?%) = 31 (! - rj) 2 {l + 277) 


3 In ^ + Yr-X r 
/ 1 2 


(1 - 77 )(34 + 2377 - 5I77 2 + I677 3 ) + 2(1 - 77) 3 (3 - 77 s ) ^ 


2(2 - 77) 

(26 — 19/7 + 4 / 7 2 )? 7 2 


(2 - >!)■ 


2 — rj 


27 


In 77 


2 — 77 


V2}. 

(A.12) 

8 < 1 - 6 ^ M. 

9 ' 

(A. 13 ) 


A. 3.2. 1 S [ q ] 


Si[ I Si 1) ]W = (- 4 )(l-tf 


81+2(77) + 4 ln(l — rj) In77 — AZr + 


47T 2 


1677(1 - 77) In77 + ^^ ln(l - 77) + 20 (1 - 77) 2 . 


5 ' 2 [ 1 *S'o 1) ]( 77 ) = 4(1 — 7 7 ) 2 


77 


2 

771 " 

31n^ + y R -A R 

/i 2 


+ 477 s In 77 


2(1 - 77 )(34 - 53/7 + I777 2 ) + 8(l- 77) 3 (3 - 77) ^ ^ 


2 — 77 


(2- v y 


gsl'S^iv) = | ((- 1 ) (1 - rj) (11 - 7?7 + 277 2 ) — 6 In77) . 


(A. 14 ) 


(A. 15 ) 
(A. 16 ) 
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A.3.3. 3 P 0 (1) 


sipPcP’lM = o. 

IftpPcf’lM = 0. 


(A.17) 
(A.18) 


SaPfflM = y (1 - ») 2 (1 + H 


— In 


9 


Ami 


+ 2 ln(l — rf) 


- (1 — 12 rj 2 + &rf) In rj - (1 — 77) (25 — 13?7 — I8/7 2 ). 


9 


9 


(A.19) 


A. 3 . 4 . 3 Pl 1] 


£i[ 3 Pi (1) ]( 77 ) = (1 - ij) 2 (1 + 2 rf) 


8 Li 2 (77) + 4 ln(l — rf) In 77 — 4 Z R + 


An 


2n 


32 16 

- — 77 (1 + rf) (1 - 2 r f) In 77 - — (1 - if) 2 (5 + Arf) ln(l - rf) 
o o 

+ ^ (1 - rf) (5 + 9 1] - 6 1] 2 ). 

o 


(A. 20 ) 


92[ 3 Pi 1) }(v) = | (1 - 9) 2 (l + 2 V ) ( 3 In ^ + Y R - X 


V 9 


R 


64 


V 


Li 2 ( 1 -+ Li 2 f --—^ - 2 Li 2 (?/) - In 2 ln (2 - rf) + In 2 In77 


2 — 77 


2 — 77 


M hl(1 _ „ + + 16(1 -’»W^ 10 ’> 3 + 3 P mi - ,) 
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(A. 21 ) 
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A.3.5. 3 P 9 (1) 
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A. 3 . 9 . 3 pj 8) 


We have summed over J = 0 , 1 , 2 . 


9 i[ 3 Pj 8 ) ](v) = (1 - h) 2 (l + 2 rj) [(- 48 ) In + 96 ln(l - rj) 
- 24 (1 - 12?7 2 + 8 rf) In rj - 4(1 - rj) (35 + 41 rj - 94 ?/ 2 ). 




(A. 35 ) 


S2[ 3 -Pf](h) = 6(1 - v) 2 (l + 2 rj) ( 31 n^ + Yr - X. 


+ 48/ 7 2 T — Li 2 +Li 2 


h 

2(1 — ??) 
2 — 77 


i? 


2 Li 2 (77) — In 2 ln(2 — 77) + In 2 In 77 


, W1 x , ^ , 12(1 — 77) (3 + 877 s — IO/7 3 + 3?7 4 ) , ^ , 

In 77 ln(l - 77 ) + y I +-7^7- y 2 - ln(l - rj) 

6 r / 2 (48 - 9677 + 59t7 2 - 12 / 7 3 ) ^ 24(1 - 77) 77 (4 - 2I77 + I777 2 - 477 s ) ^ 2 


( 2 - 77)2 

3(1 - 77) (34 + 2377 - 3t 7 2 - 877 s ) 
2 — 77 

^3[M 8) ](h) = (1 - h) 2 (1 + 2/7) (- 30 ) In 


(2 - T7) S 


(A. 36 ) 


7 T 7 2 4 

/i" 3 


+ 9 (1 - >)) 2 (1 + 2 r,) | - 4 In ^ + i.X's + yV's - | Z„ - 3 1 h 2 (2 - ,) 


2 

3 


+ 6 ln(l — 77) ln (2 — 77) ) — 12 (3 — 477) (3 + 7 r ]) ( Li 2 ( — In 2 ln (2 — 77) 


+ In 2 In 77) + 6 (29 + 3677 — 9I77 2 — I477 3 ) Li 2 (77) 


+ 6 (9 + I8/7 — 29 ? 7 2 — I8/7 3 ) Li 2 


1 — 77' 

2 — 77, 


- (7 + 3677 - 25 t 7 2 - 5877 s ) vr 2 


+ 12 (5 + 977 — I677 2 — 877 s ) In 77 ln(l — rj) 

30(1 - 77 )(14 + 677 - 69 t 7 2 + 58 t 7 3 - I377 4 ) 

(2 — rj ) 2 

3(4 - 6877 - 15977 2 + 632 ? 7 3 - 466 t 7 4 + IO6/7 5 ) 


(2 - 77 y 


ln(l — 77) 
In 77 


33 



(A.37) 


6(1 - 77)(36 + II277 - 399/7 2 + 269t7 3 - 58t/ 4 ) 

-7-77^- In 2 

(2-77)2 


(1 - 77 ) (1078 + 85177 - 3757 T 7 2 + 140 2 / 7 3 ) 

+ W^) 


A.3.10. 1 P 1 (8) 


5'i[ 1 ^ > i 8 ) ](^7) = 16 (1 — 77) s 


— In 


b 

4 m 2 


+ 2 ln(l — 77 ) 


- (1 — 677 + I 277 2 ) In 77 -(1 — 77 ) (119 — 8577 + 8rf). 

3 9 


92[ 1 Pi 8) ](v) = 0 


33 M * 1 ] (' 7 ) = 10(1 - 77) 5 


In 


jr_ 

4t77 2 


+ 2 ln(l — 77 ) 


- - (7 - 1577 + 30?7 2 ) In 77 - - (1 ~ 77 ) (347 - 244t7 + 29t7 2 ) 
3 9 


(A.38) 
(A.39) 


(A.40) 


APPENDIX B. SHORT-DISTANCE COEFFICIENTS FOR ENERGY 

DISTRIBUTIONS 

At leading order in a s the energy of the quarkonium is fixed. A non-trivial energy 
distribution is generated by gluon emission at NLO. The functions gi[n](rj, x), where x 
is the quarkonium energy fraction as defined in the text, can be expressed in the form 

gi[n](v,x) = gi[n](rj) 5(1 + 77 - x) + \gi, rea i[n\(ri,x)\ + , (B.l) 

where g l [n] ( 77 ) is the corresponding function for the total integrated decay rate, given in 
Appendix A.3. The ‘plus-distribution’ is defined by 

1+»J l+»7 

J dx[f(x)] + t(x) = J dxf(x)(t(x)r*t(l+ri)) (B. 2 ) 

2 vA 2 Vv 

for a test function t. We also introduce 

p = \Jx 2 — 477 . (B.3) 

The kinematic limits on x are < x < 1 + 77 . In the following we give the energy 
distribution functions (ji^eai[n](:V, x )- 
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g 3 ,real[ 3 Si ] ](x, 7 ]) = ^ (2 - x)p + y x(l + V - x) 111 ^ + ^ • 


(B.4) 


(B.5) 


(B.6) 


B.2. 's,\ 


fi'i,rea/[ 1 5 'o } ](a;, 77 ) = - 


28(1 - 77) 

l + ?7 — a; 


- 12 p 


I + 77 — 2; / 2 — 2; — p 


fi'2,reaz[ 1 *S'J 1) ](^, '7) = ~8 p + 8(1 + 77 -x) 111 


2 — 2; + p 

2 — 2; — p 


fi'3,reai[ 1 *S' ( S 1) ](2;, 'l/) = (-4) (6 + 8/7 - 72 ?) p + 16 (1 + 77 - x) (x - 277) In ^. 

2 77 — 2; + p 


(B.7) 


(B.8) 


(B.9) 


B.3. 3 1+ 


tfl.reaiP^K 3 ^) = 0. 
^.realP-Fo 15 ]^^) = °- 


(B.10) 

(B.ll) 


g3,real( 3 Pi 1] }(x, 7/) = ^ 1 ^ ^ (26 - 40?/ + llx) j p 

+- ( 2r/ 2 + x - 772; - x 2 ) In 2r l- x -P ' 

3 v 7 2?/ — x + p 


(B.12) 
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B. 4 . ' T , 


3 r>(l) 


g^real^Pl^KXiV) 


8 \ /7(1 — 77 X 1 + 2 / 7 ) 


3 J \ 1 + 77 — x 

16 / 2(1 — 77) 2 (1 + 277 ) 


+ 3 - 2 77 p 


3 \ 1 + 77 — x 

r3 D (i)i^ ^ _ 16 f 4r > 2 


+ 3 + 577 — 2?f — x — 2rjx In 


2 — x + p 
2 — x — p 


(B.13) 


92,real [ P{ } (x, Tj) = y J + "'_ ~ “ (1 “ J P 

16 ( Arf \ 2 - x + p 

+ T (iTi+x + 1 ~~ 3 _ ) 111 2^7 


32 


3 \ 1 + 77 — x 


2 " 2 + 37)l n 2, '~ I ~7 

J 2rj — x + p 


93,real[ 3 Pl 1] }(x , Tj) = jj “ ( 16 + 34 P “ 23a; )^ P 


— — ( 3/7 + 7n 2 — 2x — 6nx + 2a: 2 ) In —— 
3 V ' 2 n — 


2r) — x — p 


2rj — x + p 


(B.14) 


(B.15) 


B. 5 . < ¥t 


3 p(l) 


= 0- 

92,real[ 3 P2 1) }{x,ri) = 0 . 

P 3 ,rea/[ 3 -P 2 (1) ](^, v) = ^ - (4 - 277 + x) j p 

+ — (377 — 5 ?7 2 + 2 a: — 277a: — 2 a: 2 ) In ——-- 

15 v 7 2 /7 — a: + p 


(B.16) 

(B.17) 


(B.18) 


B.6. IP, 


ip(4 


= 0. 

g2,real[ 1 Pi 1) ](x,Tl) = 0. 


(B.19) 

(B.20) 
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g^real^Pl^Xi Tl) = ^ ^ ~ (14 + p 

- — (2 77 -x- V x + x 2 ) In 2>1 - - - 9 . (B.21) 

3 2 1) — x + p 

B.7. 3 s{ 8) 

gi,reai[ 3 Si 8} ](x, r)) = 4(2 - x) p + 16x(l + rj - x) ln^ ——— (B.22) 

2rj — x + p 

g2,real[ 3 S[ 8) }{x, 7]) = ~2 (1 - 2 7j) p + 2 (l - 3^ - x) 111 | — ' 7 ^ - 

z oc p 

+ 8, 2 In — ~ x ~ P (B.23) 

' 2t j-x+p ' 

93,™jpSi 8) ](i,,) = ^- 32(1 ~ + ’lKl + a ? ) + 39 + 68, - ^ p (B.24) 

4 - ( - 2(1 - ’ , )( 1 + 2, ')( 5 + 4 ’') + 17 + 5 , - 8, 2 - 9 * - 8 , 7 ] In 2 -* + £ 

\ 1 + Tj — X J 2 — x — p 

+ ( 18(1 - ,,)(1 + 2,,) - 18 + 18, + 92, 2 + 10* - 26,* - lO* 2 ) In 2 ’'-* ~ P . 
y 1 + rj — x ) 2 r\ — x + p 


B.8. '5q 8) 


fl'i,reai[ 1 >S'(S 8) ](a:, 77) = -12 (6 + 877 - 7 a) p - 48(1 + 77 - x)(2// - x) In ---— 

2r) — x + p 

g2,real[ 1 So ) ](x, 7/) = -6/0 + 6(1 + Tj ~ x) 111 \. 

2 — x — p 


g 3 ,reai[ 1 S^ ) ](x,ri) = - 51 - 6077 + P 

( 6(1 -77) (5 + 477) 


1 + 77 — X 
/ 54(1 — 77) 


+ 51 + 377 — 27a; In 


2 — x + p 
2 — x — p 


\ 


1 + 77 — a; 


54 — 677 — 6O/7 + 30a; + 9077a; — 30a; In 


2 rj — x — p 
2 r ] — x + p 


(B.25) 

(B.26) 

(B.27) 
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We have summed over J = 0,1,2. 


B.9. 6 P 


3 d( 8 ) 


gi,reai[ 3 Pj b) ](x,ri) = l^ 96 ^ - 12(18 + 16// - 13a;) j p 

+ 48 (—2// — 8// 2 + 3a; + 5 px — 3a; 2 ) In —————. 
v 7 2// — x + p 


(B.28) 




(i^r 12 ^ 2 ^ 


( 48r/ 


\ 1 + T) — X 

48// 2 

1 + 77 — a; 


2 +12(l-3„-U) l, d- X + P 

J 2 — x — p 

2 // — x — p 


+ 72 // 2 In 


2p — x + p 


(B.29) 


93,«<j[ 3 d 8) ](i,>?) = (~ 12(11 i+l V _ X 32V * ) - 24 < 5 - 12 '< + *)) P (B- 3 °) 

+ f- 12(5 + 9, '- 16, ' 2 - 8 ' ,3) + 6(17 - 4,, - 8„ 2 - 9* - 8„7) In 2 ~ X + P 
\ I + 77 — a; J 2 — x — p 

+ ( 12(3 ~ 4 " )(3 + 7,,) - 6(18 - 8 , - n,f - 24x + 2„x + 24* 2 )) ,n 

y 1 + 77 — x J 2// — x + p 


B.10. 1 P 1 (8) 

^rea/PA^K^) = ~ 4 ( X4 + j P 

— 16 ( 2 n — x — nx + x 2 )ln— -———. (B. 31 ) 

v 7 2// — a; + p 

g 2 ,real[ 1 Pi 8 ) ](x,r J ) = 0 . (B. 32 ) 

S , 3,rea/[ 1 A (8) ](^, h) = _ >? j, ~ (62 + 36r/ - 29 a;) j p 

— 4 (l4r/ + 9// 2 — 7 a; — 16r/a; + 7 a; 2 ) In —— -— . (B. 33 ) 

v 7 2 // — x + p 
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